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The present work is based on the nonequilibrium perturbative formalism. There the self-energies are derived 
up to the forth-order. In consequence, it proves that the nonequilibrium ( real-time ) perturbative expansion 
can be connected with the Matsubara imaginary-time perturbative expansion for equilibrium. As the numerical 
results, the Kondo resonance still disappears for bias voltage exceeding the Kondo temperatures, as observed 
in experiments of two terminal systems. 
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1. Introduction 

1.1. Nonequilibrium Perturbative Formalism 

The basic idea on the nonequilibrium perturbation theory grounded on the time-contour which 
starts and ends at t — — oo via t — oo has been proposed by Schwinger. [|T] After that, the frame 
of the nonequilibrium perturbation theory has been built up using the nonequilibrium Green's 
functions given after the time-contour by Keldysh. [|2] The perturbative equation is expressed in 
matrix form: 



G = 



g 



g S G, 



(1) 



where 



G = 



G— 
G> 



G< 
G ++ 



s = 



s— 



The nonequilibrium Green's functions are given in the Heisenberg representation by 

G—{t x ,t 2 ) = -i{Td{t x )${t 2 )), 
G ++ (t u t 2 ) = -iifdihJcVih)), 
G > (t 1 ,t 2 ) = (i 2 )), 
G < (ti,t 2 ) = i($(t 2 )d(h)). 



(2) 
(3) 
(4) 
(5) 



Here, the time ordering operator T arranges in chronological order and T is the anti time ordering 
operator which arranges in the reverse of chronological order. The angular brackets denote thermal 
average in nonequilibrium. 

The present work is undertaken on the basis of the nonequilibrium perturbative formalism, Eq. 
( 1 ). The retarded and advanced self-energies up to the fourth-order are formulated. Then it is 
confirmed that the nonequilibrium perturbative expansion can be connected with the Matsubara 
imaginary-time perturbative expansion for equilibrium. 
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1.2. The Kondo effect 

The Kondo effect [ [3] was discovered forty years ago and after that, the Kondo physics has 
been clarified from Landau's Fermi liquid theory [ |4] , the renormalization group [ [5] , scaling [ [6] , 
etc.. Besides, generalized Kondo problem, that of more than one channel or one impurity has been 
investigated. [El [8] Then, the Kondo effect in electron transport through a quantum dot has been 
predicted theoretically at the end of 1980s [ [9] and after a decade, this phenomenon has been 
observed. [|TU] The Kondo effect has been studied theoretically using the Anderson model and the 
predictions have been confirmed experimentally. In the Kondo regime, the conductance has been 
observed to reach the unitarity limit and the Kondo temperatures estimated from observation [ lllj 
are in excellent agreement with the expression derived using the Anderson model. [ 112] Furthermore, 
the Kondo effect in a quantum dot has been studied for nonequilibrium system where the bias 
voltage is applied. [ [T3] The Yamada-Yosida theory, [ [T3] perturbation theory for equilibrium 
based on the Fermi liquid theory [[4] has been extended to nonequilibrium system and it has been 
shown that for bias voltage higher than the Kondo temperatures, the Kondo resonance disappears 
in the spectral function with the second-order self-energy of the Anderson model. [ [15] This is in 
good agreement with the experiments of two terminal systems that it has been observed that the 
Kondo effect is suppressed when source-drain bias voltage is comparable to or exceeds the Kondo 
temperatures. [ [TBI [T7] 

In the present work, using self-energies derived up to the fourth-order, the behavior of Kondo 
resonance is investigated for nonequilibrium state caused by bias voltage. The numerical results 
are still that the Kondo resonance is broken, as supported by two-terminal experiments. 



2. Nonequilibrium Perturbation Theory 
2.1. Formalism 

A thermal average can be obtained on the basis of the nonequilibrium perturbation theory. [ 
[H [21 [TBI US EQl Ell [22] It is assumed that we can know only the state at t — — oo, that is, initially 
at t = — oo, the system is equilibrium and/or noninteracting state. The perturbation is turned on 
at t = — oo and introduced adiabatically and then, brought wholly into the system at t = 0; around 
t = 0, the system is regarded as stationary nonequilibrium and/or interacting state. After that, 
the perturbation is taken away adiabatically and vanishes at t — oo. When the time evolution of 
the state is irreversible, then, the state at t — oo cannot be well-defined. The time evolution is 
therefore, performed along the real-time contour which starts and ends at t = — oo, as illustrated 
in Fig. 1. 



+branch 
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OO 



-OO 
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Figure 1. The time-contour which starts and ends at t = — oo. 
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S matrix is defined by 
SO, to) = 



! V h 



dt-L, 



dt„T 



n l (t 1 )...n 1 (t n ) 



= T 



cxp <j — / dt Hi{t ) 



5(t,t ) t - 5(i ,t)=f 



exp 



dt Hi(t 



(G) 
(7) 



Here Tli is perturbation term in interaction representation. 

For thermal equilibrium, the statistical operator ( density matrix ) is written in Gibbs form for 
the grand canonical ensemble by 



QG 



Tre-/3(W-MiV) 



(8) 



Equation ( 8 ) is not valid exactly for nonequilibrium. We have no specific limitations upon the sta- 
tistical operator. The statistical operator can generally be expressed in Schrodinger representation 
by[[2Ql[22] 



Qs® = E l TOs ^ > Pm < m s(*)l- 



(9) 



Here, P m is probability that the system is in state m and \ms(t) > is the state in Schrodinger 
representation. g$ satisfies the Liouville equation by 



(10) 



The statistical operator in the interaction representation is given by g(t) — e lUat l h Qs(t)e in ot/h 
and satisfies 



ih 



dp 

'at 



[Hi,q]. 



(11) 



As a matter of course, qs(0) = g(0) = g(0). Here g(t) is in the Heisenberg representation. The 
time evolution is expressed using S matrix by 



g(t) = S(t,t )g(t )S(t ,t). 



(12) 



The thermal average in the Heisenberg representation at t — can be obtained, for example 
by [HUES] 

(TA(t)B(t )) 
= Tr[g(0)T A(t)B(t')} 
= Tr[e(-oo)5(-oo,0)TA(t)5(i')5(0,-oo)] 
= Tr[e(-oo)5(-oo,oo){T5(oo,-£»)i(t _ )B(i'-)}] 



oo oo 

EE 

n— 1 m— 1 

X 



1 1 

n\m\ \h 



dt 1 . 



dt. 



— oo J —oo 



dU 



dt, 



{f W x (t+). . {THxCtr). . .Hi(t-)A(t-)B(t'-)} ) av , 

where (•••)„-,= Tr[g(— oo). . .]. Then, the thermal average is derived by following the ordinary 
procedure via the Wick's theorem. 
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2.2. Relation of Self-Energy 



After the perturbative expansion is executed, the retarded and advanced self-energies are for- 
mulated. According to the definition, the retarded and advanced Green's functions are given by 



G r {hM) = - iO(ti - t 2 )({d(t 1 ), S{t 2 )}), 
G a (h,t 2 ) = iO(t 2 - *i)<{d(ii), d\t 2 )}). 



(13) 
(14) 



Here, the curly brackets denote anticommutator. The Dyson's equations for the retarded and 
advanced Green's functions are given by 



G" 
G" 



!J 



g r £ r G r , 
g a £ a G a . 



(15) 
(16) 



As the necessity to Eqs. ( 15 ) and ( 16 ), the self-energies £ r and E a are also required to be 
retarded and advanced functions in time, respectively. In accordance with the ordinary procedure 
of nonequilibrium perturbative formalism, [|2| I19j there 



1 



1 -1 
1 1 



and using this, then, 



£ = 



s— 

E> 



E< 



LSLt 



Q £ r 
£ Q 



The relationship for self-energies ought to be obtained here by comparison of Eq. ( 1 ) with Eqs. 
( 15 ) and ( 16 ): 



17 (t) 
E Q (t) 



£" 



•(t) + E<(t) 
"(t) + E>(t) 



-£ 
-E H 



"(*)-£>(*), 
-(i)-S<(t). 



(17) 
(18) 



3. Expressions of Self-Energy for Anderson model 
3.1. Anderson Model 

We consider equilibrium and nonequilibrium stationary states. Nonequilibrium state is caused 
by finite bias voltage, that is, the difference of chemical potentials; after bias voltage was turned on, 
long time has passed enough to reach stationary states. Since the states are stationary, Hamiltonian 
has no time dependence. The system is described by the Anderson model connected to leads. The 
impurity with on-site energy Eq and the Coulomb interaction U is connected to the left and right 
leads by the mixing matrix elements, vl and vr. The Anderson Hamiltonian is given by 

H= E ^2 nda + ^ ™Lo + Mfl fi Ra + U (n d} - (n d f)) (n di - (n dl )) 

a a a 

-J2vL(dlc La +H.c.) -^ R (4c HCT +H.a). (19) 

& (d) is creation (annihilation) operator for electron on the impurity, and c L and (c^ and c^j) 
are creation (annihilation) operators in the left and right leads, respectively, a is index for spin. 
The chemical potentials in the isolated left and right leads are y,L and /i/j, respectively. The applied 
voltage is, therefore defined by eV = fii — M_r- 

We consider that the band-width of left and right leads is large infinitely, so that the coupling 
functions, Tl and Tr can be taken to be independent of energy, E. On-site energy Eq is set cancel- 
ing with the Hartree term, i.e. the first-order contribution to self-energy for electron correlation: 
K il \E)=K {1 \E)=U(n_ a ). 
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Accordingly, the Fourier components of the noninteracting ( unperturbed ) Green's functions 
reduce to 

f(E) = ^_ (20) 

9«(E) = ^ (21) 

where T = (Tl + r^)/2. Hence, the inverse Fourier components can be written by 

g r (t) = -i6(t)e~ T \ 
g a {t) = i6(-t)e rt . 

In addition, from Eq. ( 1 ), we have 

9 <{E) = g r (E) [if L (E)T L +if R (E)T R ] g a (E), (22) 
g>(E) = g r (E) [i(f L (E) - 1)T L + i(f R (E) - 1)T R ] g a {E). (23) 

/i and f R are the Fermi distribution functions in the isolated left and right leads, respectively. By 
Eqs. ( 22 ) and ( 23 ), the nonequilibrium state is introduced as the superposition of the left and 
right leads. Then, the effective Fermi distribution function can be expressed by [ 115] 

i-L+i-R 

3.2. Self-Energy 

The retarded and advanced self-energies are derived up to the forth-order. Equations ( 17 ) and 
( 18 ) are divided into retarded and advanced terms in time: 

£ r (i) = [£— (t)+£<(t)]0(t) + [E— (t)+E<(i)]0(-t) 

= -[£++(«) + E>(t)]0(t) - [£ ++ (t) + E>(t)]0(-i), 

£ a (i) = [E— (t) + £>(t)]0(t) + [£— (t)+£>(t)]0(-t) 

= -[£++(«) + E<(t)]0(t) - [E++(t) + E<(t)]f(-t). 

Then it is found that for the self-energies drawn using the Wick's theorem, 

E (t)9(t) = -£>(*)#(*), £++(t)0(t) = -E<(t)0(i), 
E (t)0(—t) = -£<(t)0(-t), E++{t)6(-t) = -£>(i)0(-t). 

It leads to 

[E (<) + £<(t)]0(-t) = -[£++(*) + E>(t)]0(-t) = 0, 
[E (i) + E>(t)]0(t) = -[£++(*) + £<(*)]#(*) = 0. 

As a consequence, the retarded and advanced self-energies are obtained as retarded and advanced 
functions of time, respectively: 

£ r (i) = (t) + E<(t)]0(t) = -[E ++ + (t)]0(t), 

£ a (i) = [£— (t) + £>(f)]0(-t) =-[£++(*) +£<(*)]#(-<). 



In addition, it proves 



£'<-(*) = [£<(<)- £>(i)]0(t), 
£ Q (t) = [£>(*) -E<(t)]0(-f). 
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Hence the following relations still stand: £ r — £ a = £ 

G< = (l + G r £ , ') 5 < (l + G a £ a )-G r £ < G a : 
G> = (1 + G''£ r )5 > (l + G Q £ a ) - G r £>G a 



£> and furthermore 



As the results, the second-order self-energy is written by 



DO 



2 / dhe* Etl 
'o 



dt\e 



iEt-i 



g>(h)g>(t 1 )g<(-t 1 ) 

9 ± (ti)9 > (ti)g<(-t 1 ) 
+9 < (h)g ± (ti)g > (-t 1 ) 
+9 < (h)g > (ti)g ± (-t 1 ) 



(25) 



£«(2)(.E) = 



iEli 



dt\e 



iEt x 



9 <(t 1 )g<(t 1 )g>(-t 1 ) 

_ -9 > {h)g > {tx)g < {-tx) 

5 ± (*i)ff > (ti)5 < (-*i) 
+3 < (ii)ff ± (^i)ff > (-ii) 
+3 < (ii)3 > (^i)3 ± (-ii) , 

VI 



(26) 



Here £/ ± (t) = g r (t) + g a {t), that is, g+{t) = g r {t) = -i0(i)e~ r * for t>0 and = g a (t) = 

id(-t)e rt for t < 0. Additionally, g<(t) and g>(t) are the inverse Fourier components of Eqs. ( 22 
) and ( 23 ). Figure 2 shows the diagram for the second-order self-energy. As shown numerically 
later, the second-order contribution coincide with those derived by Hershfield et al. [HI]. In the 
symmetric equilibrium case, the asymptotic behavior at low energy is expressed by 



the exact results based on the Bethe ansatz method. [23, 24 



E 



II 



E 



(27) 




Figure 2. The diagram for the second-order self-energy. The solid line denotes the noninteracting 
Green's function and the dashed line denotes interaction. 

The third-order terms corresponding to the diagram in Fig. 3(a) are expressed by 



Kf\E) = 



I dti 
'o 



dt 



2C 



iEt x 



g < (-t 1 )g>(t 1 -t 2 )g>{t 1 -t 2 ) 
-g > (-t 1 )g<{t 1 -t 2 )g<{t 1 -t 2 ) 



x g ± {t 2 )g>{t 2 )+g<{t 2 )g ± (t 2 ) , 



(28) 



dt 2 e iEtl 



g>(-t 1 )g<(t 1 -t 2 )g<{t 1 -t 2 ) 
~g < (-t 1 )g > {t 1 ~t 2 )g > (t 1 -t 2 ) 

x [ g ± (t 2 )g>(t 2 )+g<(t 2 )g±(t 2 ) ]. 



(29) 
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Figure 3(b) illustrates the diagram for the following terms: 



^(3) 
-'ph 



Iq 



dt%e 



iEti 



g>(t 1 )g>(t 1 -t 2 )g<(t 2 -t 1 ) 
-g < {t 1 )g<(t 1 -t 2 )g>(t 2 -t 1 ) 



x [ g ± (t 2 )g<(-t 2 ) + g<(t 2 )g±(-t 2 ) 



-,o(3) 



(E) = U 



dti 



dt 2 e 



iEti 



g < (h)g < (t 1 -t 2 )g>(t 2 -t 1 ) 

-9 > {h)g > {tx-h)g < (h-t 1 ) 
x [ g ± (t 2 )g<(-t 2 ) + g<(t 2 )g±(-t 2 ) ]. 



(30) 



(31) 



Equations ( 28 )-( 31 ) for equilibrium agree exactly with those derived from the Matsubara 
imaginary-time perturbative expansion for equilibrium and analytical continuity by Zlatic et al. [ 
125] As mentioned later, it is numerically confirmed that the third-order contribution vanishes for 
the symmetric Anderson model; this is in good agreement with both the results derived from the 
Yamada-Yosida theory [ [T31 [2H and obtained on the basis of the Bethe ansatz method. [ [53] 





Figure 3. The two diagrams for the third-order self-energy. Left:(a) and Right:(b) 



ooo C?, (O OyO O fP> 



CS> O O 



O 



M O O 



Figure 4. The twelve terms for the proper fourth-order self-energy divided into four groups: 
(a)-(c), (d)-(f), (g)-(i), and (j)-(l). 



Furthermore, the fourth-order contribution to the self-energy is formulated. ( See Appendix. 
) The twelve terms for the proper fourth-order self-energy can be divided into four groups, each 
of which comprises three terms. The four groups correspond to the diagrams denoted in Figs. 4 
(a)-(c), Figs. 4 (d)-(f), Figs. 4 (g)-(i), and Figs. 4 (j)-(l), respectively. For symmetric Anderson 
model at equilibrium, the asymptotic behavior at low energy is approximately in agreement with 
those based on the Bethe ansatz method f [23] : 



(32) 
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4. Numerical Results and Discussion 

4.1. Self-Energy 

The third-order terms, Eqs. ( 28 )-( 31 ) cancel under electron-hole symmetry not only at 
equilibrium but also at nonequilibrium: S^ 3 '(E) = —J^\E) and Yf p ^(E) — ~T,p P 3 \E). As a 
consequence, the third-order contribution to self-energy vanishes in the symmetric case. It agrees 
with the results of Refs. [14, 24, 26] based on the Yamada-Yosida theory that all odd-order contri- 
butions except the Hartree term vanish for equilibrium in the symmetric single-impurity Anderson 
model; probably, it is just the same with nonequilibrium state. On the other hand, the third-order 
terms contribute to the asymmetric system where electron-hole symmetry breaks and furthermore, 
the third-order terms for spin-up and for spin-down contribute respectively when the spin degen- 
eracy is lifted for example, by magnetic field. For the fourth-order contribution, three terms which 
constitute each of four groups contribute equivalently under electron-hole symmetry. Moreover, 
to the asymmetric system, the terms brought by the diagrams of Figs. 4(a) and 4(b) contribute 
equivalently and the terms by the diagrams of Figs. 4(j) and 4(k) make equivalent contribution, 
and the rest, the eight terms contribute respectively. Further, the twenty-four terms for spin-up 
and spin-down take effect severally in the presence of magnetic field. 

The second-order and the fourth-order contributions to self-energy for zero temperature sym- 
metric Anderson model are shown in Figs. 5(a) and 5(b) and in Figs. 6(a) and 6(b), respectively. 
Equation ( 27 ) represents the curves around E = denoted by solid line in Figs. 5(a) and 5(b), 
respectively, and Equation ( 32 ) represents approximately those shown in Figs. 6(a) and 6(b), 
respectively. The curves of the second-order self-energy shown in Figs. 5(a) and 5(b) are identical 
with those of expressions derived by Hershfield et al. [HI]. In comparison of Figs. 6(a) and 6(b) 
with Figs. 5(a) and 5(b), it is found that the fourth-order contribution for equilibrium has the 
same but narrow curves at low energy with those of the second-order contribution. In addition, 
the broad curves are attached at high energy for the fourth-order self-energy. ( The higher-order 
contribution is, the more the curves must oscillate as a function of energy. ) When the voltage, 
eV/T exceeds ^2.0, the behavior of curves of self-energy changes distinctly and comes to present 
striking contrasts to that for the second-order contribution. Especially, the curve for the imaginary 
part of the fourth-order contribution rises up with maximum at E — 0. On the other hand, for 
the second-order contribution, a valley appears with minimum at energy of zero— it is quite the 
contrary. Moreover, from these results, it is expected that the sixth-order contribution to imaginary 
part of self-energy has minimum at E = 0. Because of these, the perturbative expansion is hard to 
converge for eV/T > ^2.0, as mentioned later. 

Besides, the current conservation is mentioned. In Ref. [15], it is shown that the continu- 
ity of current entering and leaving the impurity stands exactly at any strength of U within the 
approximation up to the second-order for the symmetric single-impurity Anderson model. In com- 
parison of Figs. 6(a) and 6(b) with Figs. 5(a) and 5(b), it is found that curves of fourth-order 
self-energy have the symmetry similar to those of the second-order. From this, it is anticipated 
that the current conservation are satisfied perfectly with approximation up to the fourth-order in 
the single-impurity system where electron- hole symmetry holds. The continuity of current can be 
maintained perfectly in single-impurity system as far as electron-hole symmetry stands. On the 
other hand, current comes to fail to be conserved with increasing U in asymmetric single-impurity 
case and in two-impurity case. 

4.2. Spectral Function 

The spectral function with the second-order self-energy is generally known. It is plotted for 
U/T = 10.0 and zero temperature in Fig. 7. For equilibrium, the Kondo peak at energy of zero is 
very sharp and the two-side broad peaks appear at E~ ±{7/2. The curve for eV = is identical 
with that shown in Ref. [24]. As eV becomes higher than the Kondo temperatures, ksTx [ I27j. 
the Kondo peak becomes lower and finally vanishes, while the two-side broad peaks rise at E~ 

±U/2. [ng 
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Figure 5. The second-order self-energy for the symmetric Anderson model at U/T — 1.0 and 
zero temperature, (a) The real part and (b) the imaginary part at equilibrium ( solid line ), 
eV/T = 1.0 ( thin solid line ), and eV/T = 2.0 ( dashed line ). 




Figure 6. The fourth-order self-energy for the symmetric Anderson model at U/T = 1.0 and 
zero temperature, (a) The real part and (b) the imaginary part at equilibrium ( solid line ), 
eV/T = 1.0 ( thin solid line ), and eV/T = 2.0 ( dashed line ). The fourth-order contribution 
for equilibrium has the same but narrow curves at low energy with those of the second-order 
contribution. 
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0.8- 

g 0.6- 
O 

g 0.4- 




E/r 

Figure 7. The spectral function with the second-order self-energy at U/Y = 10.0 for the symmet- 
ric Anderson model at equilibrium ( solid line ), eV/Y = 1.0 ( thin solid line ) and eV/Y — 2.0 
( dashed line ). 




E/T 

Figure 8. The spectral function with self-energy up to the fourth-order at equilibrium for the 
symmetric Anderson model at U/Y — 3.5 ( dashed line ) and U/Y — 5.0 ( solid line ). 
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Figure 8 shows the spectral function with the self-energy up to the fourth-order for equilibrium 
and zero temperature. With strengthening U, two-side narrow peaks come to occur in the vicinity 
of E = ±U/2 in addition to the Kondo peak. At U large enough, the Kondo peak becomes very 
acute and two-side narrow peaks rise higher and sharpen; the energy levels for the atomic limit 
are produced distinctly. The fourth-order self-energy has the same but narrow curves as functions 
of energy with those of the second-order and those curves make the peaks at E = ±U/2. 

For the present approximation up to the fourth-order, the Kondo peak at E = reaches the 
unitarity limit and the charge, (n) corresponds to 1/2, that is, the Friedel sum rule is correctly 
satisfied: [[28] 

p(E f ) = sin 2 (^(n))/7rr, (33) 

where p{Ef) is the local density of states at the Fermi energy. Here, the discussions should be 
made on the ranges of U in which the present approximation up to the fourth-order stands. From 
the results, it is found that the approximation within the fourth-order holds up to U/T ~5.0 and 
is beyond the validity for U/T> ~6.0. In addition, the curve for imaginary part of the fourth-order 
contribution is positive partly, as shown in Fig. 6(b) and as a consequence, the curve of the spectral 
function becomes negative partly for too large U. In such a case, the present approximation is out 
of validity and the higher-order terms are required. 




E/r E/r 

Figure 9. The spectral functions with self-energy up to the fourth-order at eV/T = 0.5 ( Left 
) and eV/Y — 1.0 ( Right ) for the symmetric Anderson model. U/Y = 3.5 ( dashed line ) and 
U/T = 5.0 ( solid line ). 

Next, the results for nonequilibrium and zero temperature are shown. The expression for the 
Friedel sum rule, Eq. ( 33 ) does not stand for nonequilibrium, since the charge cannot be expressed 
with respect to the local density of states. All the same, the Kondo peak reaches the unitarity limit 
and (n) = 1/2 in the symmetric and noninteracting case. The spectral functions with the self- 
energy up to the fourth-order are plotted for eV/T — 0.5 and eV/T = 1.0 in Figs. 9, respectively. 
When U is strengthened and eV exceeds /cbTr- ( approximately, ksTx/T ~0.5 for U/T — 3.5 
and IcbTk/T ~0.3 for U/T— 5.0 ), the Kondo peak for eV/T = 0.5 falls in and instead, the two- 
side narrow peaks remain to sharpen in the vicinity of E = ±17/2. For eV/T = 1.0, the Kondo 
peak becomes broad and disappears for U large enough. The two-side peaks is generated small in 
the vicinity of E = ±[7/2. The Kondo resonance is quite broken for bias voltage exceeding the 
Kondo temperatures; this accords with the experimental results of two terminal systems that the 
Kondo effect is suppressed when source-drain bias voltage is comparable to or exceeds the Kondo 
temperatures, eV>kBTK- [HUE] For eV/T > ~2.0, the Kondo peak does not lower even when 
eV is much larger than UbTk- The perturbative expansion is hard to converge on account of the 
imaginary part of the self-energy for eV/T > ^2.0, as described before; thereby, the higher-order 
contribution to self-energy is probably needed for high voltage. 
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In the present work, nonequilibrium state is represented as the superposition of the two leads 
and the effective Fermi distribution function, Eq. ( 24 ) is qualitatively similar to that for finite 
temperatures. From the analogy in the Fermi distribution function, it is inferred that there are 
nonequilibrium fluctuations similar to thermal fluctuations. [ I29j Because of the effective Fermi 
distribution function, not only for the second-order but also for the fourth-order, the Kondo res- 
onance is destroyed, qualitatively the same as for finite temperatures. In contrast, if the finite 
voltage state is expressed as two localized states, the numerical results of the Kondo peak splitting 
can be obtained. All the same, for two terminal systems, the Kondo resonance splitting may not 
take place for finite bias voltage. 

Summary: The present work is based on the nonequilibrium perturbative formalism. Here the 
self-energies are derived and then it is indicated that the nonequilibrium ( real-time ) perturbative 
expansion can be related to the Matsubara imaginary-time perturbative expansion for equilib- 
rium. As the numerical results, the Kondo peak disappears as bias voltage exceeding the Kondo 
temperatures. Because of the analogy of the effective Fermi distribution function for nonequilib- 
rium with that for finite temperatures, the present result is qualitatively similar to that for finite 
temperatures. This characteristic appears in the experiments of two terminal systems. 
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Appendix 

The twelve terms for the fourth-order contribution can be divided into four groups, each of 
which is composed of three terms. The four groups are brought from diagrams denoted in Figs. 4 
(a)-(c), Figs. 4 (d)-(f), Figs. 4 (g)-(i), and Figs. 4 (j)-(l), respectively. The terms for the diagrams 
illustrated in Figs. 4(a) and 4(b) are equivalent except for the spin indices and expressed by 



r(4) 



(E) 



U 4 



dti I dt2 

o J -oo 

g < {h)g<{t 1 -t 2 
-9 > (h)g > (ti - 1 2 - t 3 )g<{-h + *2 + 1 3 ) 
g ± (t 2 )g < (-t 2 ) + g<(t 2 )g ± (-t 2 ) ] 

g ± (t 3 )g<(-t 3 )+g<(t 3 )g±(-t 3 ) 1, 



dt 3 e lEtl 

h)g > {-ti+t 2 + U) 



(34) 



„o(4) 



(E) = U 4 



dh 



dte 



g > {t 1 )g>(t 1 -t 2 



dt 3 e lEtl 

X) 

t 3 )g < (-ti + t 2 + t 3 ) 



-g < {t l )g<{t l -t 2 - t 3 )g>(-t 1 + t 2 + t 3 ) 



x [ g ± (t 2 )g<(-t 2 )+g<(t 2 )g±(-t 2 ) 
x [ g ± (t3)g < (-t 3 )+g<(t 3 )g±(-t 3 ) 
Additionally, Figure 4(c) shows the diagram for the following terms: 



K (i) (E) = U 4 / dti / dt 2 / 

Jo J-oo J— i 



dt 



3e 



iEU 



g > (~h)g < (t 1 -t 2 - t 3 )g<{t 1 - t a - t 3 ) 
~g < {-h)g > (t l -t 2 - t 3 )g>(t 1 -t 2 - t 3 ) 

g ± (t 2 )g>(t 2 )+g<(t 2 )g±(t 2 ) ] 

g ± (t3)g > (t3) + g < (t 3 )g ± (h) ], 



(35) 



(36) 
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/0 r- oo n oo 

dh / dt 2 / dt 3 e lEtl 
-co J —oo J —oo 

9 < {~t 1 )g>{t l -t 2 - t 3 )s>(ti - t 2 - h) 
-9 > (-h)g < (h -t 2 - t 3 )g<(h - t 2 - i 3 ) 
x [ g ± (t 2 )g>(t 2 )+g<(t 2 )g±(t 2 ) } 
x [ 9 ± {h)g > {h)+g<{h)g±{h) ]. 

Next, the terms brought from diagram in Fig. 4(d) are expressed by 



r-OO r-OO r-OG 

^ T d {A) (E) = U 4 / dh / dt 2 / dt 3 , 

JO J — oo J —oo 



,iEt! 



g > (t 1 -t 3 )g>(t 1 -t 2 )g<(t 2 -t 1 ) 
-9 < (h - t 3 )g < (t 1 - t 2 )g> (t 2 - h) 

g > (-t 2 +t 3 )g>(t 3 )g<(~t 3 ) 
-g < (-t 2 +t 3 )g<(t 3 )g>(-t 3 ) 



xg ± (t 2 )sga(t 3 ) 



(37) 



(38) 



K {a \e) 



/0 y oo /*oo 

dh / <ft 2 / dt 3 e 4jEtl 
-oo J —oo J —oo 

' J < (h-t 3 )g < (t 1 -t 2 )g>(t 2 -t 1 ) 
-9 > (h - t 3 )g> (h - t 2 )g<(t 2 - h) 



x 5 ± (t 2 )sgn(t 3 ) 



9 > {-t 2 +t 3 )g>{t 3 )g<{-t 3 ) 
~g < (-t 2 +t 3 )g<{t 3 )g>(-t 3 ) 



The terms for diagram in Fig. 4(e) are written by 

/>00 /'OO />oo 

K {i) (E) = U 4 / rfti / dt 2 / <ft 3 < 

^0 J —oo J— oo 



9 > {h-t 2 )g>(t 1 -t 2 )g<{t 3 -t 1 ) 
-g < (h - t 2 )g<(h - t 2 )g>(t 3 - h) 

' 9 > (h-t 3 )g > (-t 3 )g<(t 3 ) 
-g < (t 2 -t 3 )g<(-t 3 )g>(t 3 ) 



x ff ± (t 2 )sgn(t 3 ) 



(39) 



(40) 



/0 /-OO />CO 

dh / rfi 2 / dt 3 e lEt i 
-oo J —oo J — oo 



9 < {h - t 2 )g<{h - t 2 )g>{t 3 - h) 
-g > (h-t 2 )g > (h-t 2 )g < (t 3 -h) 

' 9 > {t2~t 3 )g>(-t 3 )g<{t 3 ) 
-g < (t 2 -t 3 )g<(~t 3 )g>(t 3 ) 



x.9 ± (t 2 )sgn(t 3 ) 



In addition, Figure 4(f) denotes the diagram for the following terms: 



C/ 4 / dh I dt 2 / dt 3 e 

' J —oo J —oo 



X 



9 > (h-t 3 )g > (h-t 2 )g < (t 2 -h) 
-9 < (ii-i3)5 < (ii-i2)3 > (i 2 -ii) 

< (*3)fl < (*3)S > (t2-t3) 

■5 > (i3)5 > (i3)5 < (*2-*3) 



x 5 ± (-i 2 )sgn(t 3 ) 



(41) 



(42) 



/0 /*oo /*oc 

dti / dh / ^ 3 e 4iJtl 
-co J —oo J —oo 
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g<{h-t 3 )g<(h-t 2 )g>{t 2 -h) 
-9 > (*i - h)g> (ti - t 2 )g<(t 2 - ti) 

xa ± ( t Wt ) \ ^(^(t^ (h - t 3 ) 
xg ( i 2 )sgn(t 3 )^ ^g> ih ) g > { t 3 )g<(t 2 ~t 3 ) 

Next, the terms formulated from diagram illustrated in Fig. 4(g) are expressed by 



/■ OC pOC pOO 

X r g {4 \E) = U 4 / dh / dt 2 / dt 3 < 

J •' — oc J —oo 



,iEt! 



g>(t 1 )g>(t 1 -t 2 -t 3 )g<(t 2 -t 1 ) 
-9 < (h)g < (h-t 2 -t 3 )g>(t 2 -t 1 ) 

g > (t 2 +t 3 )g>{t 3 )g<{-t 3 ) 
-9 < {t 2 + t 3 )g<(t 3 )g>(-t 3 ) 



xg ± (-t 2 )sgn(t 3 ) 



/0 pOO pOO 

dh / dt 2 / dt 3 e lEtl 
-OO J —OO J —CO 



g<(t 1 )g<(t 1 -t 2 -t 3 )g>(t 2 -t 1 ) 
-9 > {ti)g > {ti-t 2 -t 3 )g<{t 2 -t 1 ) 

x n ± ( . ,, f x [ 3 > fe +t 3 )g>(t 3 )g<(~t 3 ) 
Figure 4(h) illustrates the diagram for the following terms: 



poo poo poo 

t^ 4 \e) = U 4 / dh / dt 2 / dt 3 , 

J J — 00 J — 00 



,iEt! 



g < {h)g < (h - h - t 3 ) g > (t 2 - h) 
-9 > (h)g > (h - 1 2 - t 3 ).g < (*2 - ti) 

" g > (t 3 )g > {t 3 )g < {-t2-t 3 ) 
-5 < (i3)5 < (*3)5 > (-i2-* 3 ) 



3 ± (t 2 )sgn(t 3 ) 



(43) 



(44) 



(45) 



(46) 



£/ 4 



dh I dt 2 I dt 3 e lEtl 



] > (h)g > (h-t2-t 3 )g < (t 2 -h) 

-9 < (h)g < (h - t 2 - t 3 )g>(t 2 - h) 

xn±(t u m (t ) \ 9 > (h)g > (t 3 ) g <(-t 2 -t 3 ) 

xg (i 2 )sgn(i 3 )^ _ g<{h)g < {h)g > { _ t2 _ h) 
Besides, the terms formulated from the diagram in Fig. 4(i) are written by 

JEti 



^ 4 \E) 



pOO p CO poo 

= U 4 dh / dt 2 / dt 3 ( 

J J —oc J —00 



9 > (-h)g < (h-t 2 -t 3 )g<(h-t 2 ) 
-9 < (-h)g > (h - t 2 - t 3 )g>(h - h) 

' 9 > {t2 + t 3 )g > (t 3 )g<(-t 3 ) 
-9 < (t 2 + t 3 )g<(t 3 )g>(-t 3 ) 



x.g ± (t 2 )sgn(t 3 ) 



zf 4 \E) 



U 4 







dh 



00 



dh 



2 UL 3 I 

— 00 J — 00 

<(—+. \ n >(+. _ _ +„^>f 



oc 



dt 3 e l£tl 



ff < (-ti)fl > (ti-t2-*3)ff > (ti-*2) 

^{-h^ih - * 2 - h)g < (h - h) 

g > (t2+t 3 )g > (t 3 )g<(-t 3 ) 
-g < (t2 + t 3 )g < (t 3 )g>(-t 3 ) 



x 5 ± (t 2 )sgn(t 3 ) 



(47) 



(48) 



(49) 
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Next, the terms for diagrams denoted in Figs. 4 (j) and 4(k) are equivalent except for the spin 
indices and written by 



/>oo />00 r-oo 

U 4 / dh / dt 2 / dt 3 e lEtl 

JO J —OO J — CO 

g>(t 1 )g<(-t 1 )g>(t 1 -t 2 -t 3 ) 
~9 < (h)g > (-iil^fti - h - i 3 ) 

9 ± (h)g > (t 3 )g < (~t 3 ) " 
xg ± (t 2 ) +9 < (t 3 )g ± (t 3 )g > (-t a ) 
+9 < (h)g > (t 3 )g ± (-t 3 ) 



(50) 



/0 oo f- oo 

- oo J —oo J —oo 



g<(t 1 )g>(-t 1 )g<(t 1 -t 2 -t 3 ) 
-g > (h)g < (-h)g > (t 1 -t 2 -t 3 ) 

g ± (t 3 )g>(t 3 )g<(-t 3 ) 
+g < (h)g ± (h)g > (-t 3 ) 
+g < (h)g > (h)g ± {-h) 



In addition, the terms for diagram illustrated in Fig. 4(1) are expressed by 



,r(4) 



f- OO f-OO pQ 

(E) = U 4 / dh / dt 2 \ 



dt 3 e lEtl 



g>(t 1 )g>(t 1 )g<(-t 1 +t 2 +t 3 ) 
-g<(t 1 )g<(t 1 )g>(-t 1 +t 2 +t 3 ) 

g ± (^t 3 )g>^t 3 )g<(t 3 ) 

+g < (-t 3 )g ± (~t 3 )g>(t 3 ) 
+g < (-t3)g > (-h)g ± (h) 



x.9 ± (-t 2 ) 



(51) 



(52) 



,a(4) 



(E) 



u 4 



dh 



dh 



dt 



3 e 



. 9 <(t 1 ).g<(t 1 ) 5 >(-i 1 + t 2 + i 3 ) 
-<7 > (*i)<7 > (*i)<7 < Hi+*2+* 3 ) , 

g ± (-t 3 )g > (-t 3 )g < (t 3 ) 
+g < (-h)g ± (-t 3 )g>(t 3 ) 
+g < (-t 3 )g>(-t 3 )g±(t 3 ) 



xg ± (-h) 



(53) 
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